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I. INTRODUCTION 



Recently, much interest has been generated to the studies of physical systems with sta- 
tistical properties which cannot be described within the concept of the Markovian random 
processes. This concerns anomalous transport in turbulent plasma, strange diffusion of mag- 
netic field lines in fusion plasma, the Brownian motion of macroparticles in complex fluids 
etc. (see, for example, ^ and references cited therein). One of the prob- 

lems challenging and most important at the same time is to describe the kinetic properties 
of such systems and to understand the role of non-Markovian effects of particle and energy 
transport in such systems. The key point of such calculations is the formulation of the con- 
sistent kinetic-type equations with regard to time-nonlocality of the collision integrals and 
the renormalization of the free-streaming particle propagator due to fluctuation influence on 



particle trajectory. In 
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11| , this problem was solved on the basis of the Green function 



method similar to that proposed in 12|. The main disadvantage of such a treatment is the 
uncertainty introduced by the specific splitting of higher correlations made on the basis of 
some physical arguments. 

The aim of the present paper is to derive the non-Markovian kinetic-type equations in a 
more consistent way on the basis of a rigorous solution of the appropriate BBGKY hierarchy 
(Bogolyubov-Born-Green-Kirkwood-Yvon hierarchy) for the generalized microscopic phase 
densities and their average values. 

The paper is organized in the following order. 

In Section II, we introduce the definitions used for the description of the evolution of the 
microscopic phase density and its generalization. In Section III, we deduce the evolution 
equations for average values of microscopic phase densities and construct a solution of initial- 
value problem of the obtained hierarchy of equations. In Section IV, with the use of the 
results obtained above, we develop a new approach to the description of the evolution of 
the average value of microscopic phase density; namely, we construct an evolution equation 
for this quantity and deduce a non-Markovian generalization of the Fokker-Planck collision 
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integral. 



II. THE EVOLUTION OF THE MARGINAL MICROSCOPIC PHASE 
DENSITIES 

We consider the system of a non-fixed (i.e. arbitrary but finite) number of identical 
particles with unit mass m = 1 in the space (in the terminology of statistical mechanics 
it is known as nonequilibrium grand canonical ensemble 13|). Every particle is characterized 
by the phase space coordinates Xi = {qi,Pi), i.e. by a position in the space G and a 
momentum G M^. A description of many-particle systems is formulated in terms of two 



sets of objects: by the sequence of observables and by the sequence of states 



A. The dual BBGKY hierarchy 



li 



In [ij] it was introduced the marginal observables known as the microscopic phase 
densities G{t) = {G^^\t), G^''\t), . . . ) of fc-ary type G('=)(t) = (O, . . . , 0, G[^\t), 
Gf\t),...). 

For example, at the initial time moment t = the sequence of marginal additive-type 
microscopic phase densities [14] has the form G^^\0) = (O, . . . , 5{^i — Xi), 0, . . . ) , where 5 is 
the Dirac (5-function. Correspondingly, the sequence of marginal observables of fc-ary type, 
A; > 1, microscopic phase densities (jlj) is given as follows: 

k k 

G(^)(0)= (0,...,0, Yl n^(^'-^*<)'0'---)- (1) 



il^...7^ifc = l 1=1 



Let Hn = X]r=i^ + J27<j=i^ili ~ 'ij) is the Hamiltonian of a n-particle system, 

j 

$(gi — Qj) is a two-body interaction potential, Y = {xi, . . . ,Xs), (a;i, . . . , V, . . . , x^) = 
(xi, . . . , Xj-i, Xj+i, . . . , Xs) =Y\xj and the brackets (■, ■) is a scalar product of vectors. 

The marginal microscopic phase densities Gi^\t) = Gi''\t, ^i, . . . ,C,k]Xi, . . . , Xg) of every 
k-cLTj type are governed by the initial- value problem of the dual BBGKY hierarchy jl^. jig 

|G«wHt(P<.|-)-i:(|-*te-,.).|->)Gfw 

~ E - ^^O' ^^^s%it, Y \ X,) (2) 



with the initial data 



|t=o=Gi'=)(0), .>fc>l. (3) 
As a case in point, we adduce the first equation of hierarchy ([2]) 

For the marginal additive-type microscopic phase density, the first two equations have the 
form 

6; x„ X,) = ( I") - E - ^^O' |-))^^'^(^' - 

Let Ar('=)(t) = (O, . . . , 0, N^\t), . . . , . . . ), > 1, is a sequence of microscopic 

phase densities of k-arj type [Sj 

n k 

where ^i, . . . , are the macroscopic variables C,i = (f rj) G M'^ x R'^ and the set of functions 
. . . > A; > 1, are the solution of the Cauchy problem of the Hamilton 

equations for n particles with the initial data Xi, . . . Then the sequence of marginal 
observables G{t) in terms of sequence (jl]) are defined by the formula jl^ 

Gi'\t, x„...,x,) = J2 ^ E ei, • • • , Y\{x,„ . . . , x,J), (5) 

where F = (xi, . . . ,Xs), s > 1. 

For example, if k = 1, i.e. in the case of an additive- type observable, we have the 
microscopic phase density [8|, 19| 



n 



Ni^Ht,^i;Xi,...,Xn) = ^6(^1- Xi{t,Xi,...,Xrr)). 



1=1 
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We note that expression (jl]) defines the one-parametric group of operators 9 t i— >■ 
Sn{t)Nn{0), i.e. 

iV('=)(t, ^1, . . . , xi, . . . , xr,) = 5„(t)iVf )(0), (6) 

where Ni^\o) is the initial microscopic phase density (jlj). 

Let Y = {xi, . . . ,Xs), X = Y \ {xj-^, . . . , Xj^_^}. For continuous functions in the capacity 
of initial data, a solution of Cauchy problem ([2])- ([3]) is defined by the following expansion 

Gi'\t, Y) = j2^_ E 2li+„(t, {Y \ X)^,X)Gi%{0, Y \ {x,,, . . . , x,J), (7) 

where the evolution operator 

Qi,Ut,{Y\X)^,X)= J2 (-1)'"-^(|P|-1)! n ^1^.1(^-^0 

P:{(Y\X)i,X}=U,X. X,CP 

is the (1 + n)th-oTdeT cumulant of groups S\Xi\{t) of operators ([6]), is the sum over all 
possible partitions P of the set {(F \ X)i,X} into |P| nonempty mutually disjoint subsets 
Xi C {{Y \ X))^,XY The set (Y \ X)i consists of one element of Y\X, i.e. the set 
Y\X = {xj^, . . . ,Xj^_^} is a connected subset of the partition P (|P| = 1). 

Then in terms of variables ^i, ■ ■ ■ ,C,k, for the marginal microscopic phase densities of fc-ary 
type G^'^^t) = (0, . . . , 0, (t), . . . , ci^^t), . . . ) we derive 

+ E / d^^M^.'^in - r.+i), ^)Gr'\t) (8) 
with the initial data 

s k 
il7^...^ifc = l i = l 

Here, 1 < r < s, and if A; = s, the marginal microscopic phase density Gi^\t) is governed by 
the Liouville equation. 

Thus, in terms of variables ^i, . . . the dual BBGKY hierarchy ([HD for marginal micro- 
scopic phase densities ([7]) is represented as the Bogolyubov equations set with respect the 
arity index k > 1, while evolution equations ([8]) have a structure of a sequence of equations 
with respect to the index of the number of particles s > k. 
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We note that, for every s, a solution of initial-value problem ©-([Q]) of equations set ([8]) 
can be represented as the following expansion: 

G'^J'\t,Ci,...,Ck;xi,...,Xs) = 

= T.^j ^^a+l---rfa+n2li+„(-t,ri,a.+l,...,a+n)G'('=+"H0), 
n=0 ■ 

where Gi''~^"'\o) is initial data ([9]), and 2ti+„(— t, Fi, ^^+1, . . . , is the (1 + n)t/i-order 
cumulant fl25l) of the groups of evolution operators fl22|) defined further. 



B. On the evolution of states of many-particle systems 

We furnish further comments about the evolution of microscopic phase densities in the 
framework of the evolution of states of many-particle systems. 

In this case, the microscopic phase densities are defined at the initial time moment A^(0) = 
(A^W(0),...,Ar('=)(0),...), where ArW(O) = (0, . . . , 0, A^f ^(0), . . . , M''^(0), . . .) and 

n k 

iVW(0)^iV('=)(0,ei,...,efc;:^i,...,^n)= W^i^i-^n)^ n>k>l. 

n^...^ifc=l 1=1 

The evolution of states is usually described by the Cauchy problem of a sequence of 
the Liouville equations for the sequence D{t) = (/, Z}i(t), . . . , . . . ) of distribution 

functions Dn(t) = Dn(t, xi, . . . , x„) 

= ( - Efe. aj> + E - *). e^))^M. (10) 

i=l 27^j=l 

^n(t)|t=o = ^„(0), n>l. (11) 

A solution of the Cauchy problem (fTO|) - (fTT]) is constructed in [13)]. Average values of the 
microscopic phase densities A^(0) are determined by the expressions 

(ArW)(t,ei,...,6) = 
1 °° 1 /" 

= (l,Z:'(0)) — / dXk+l. . .dXk+nDk+n{t,^l, . . . ,^k,Xk+l, ■ ■ ■ ,Xk+n), (12) 

n=0 ^' 

where Dk(t,^i, . . . ,^k, Xk+i, • • • , Xk+n) is the value of a solution of initial- value problem (fTOj) - 
(fTTj) at a point ^i, . . . , ^k, Xk+i, ■ ■ ■ , Xk+n- Thus, the evolution of the functions (A^*^^^)(t), k > 



1, will be governed by the BBGKY hierarchy with respect to the arity index k (see next 
section) . 

In the thermodynamic limit, the evolution of states is described in terms of the marginal 

nn 

distribution functions governed by the BBGKY hierarchy 113|,[16[. From Liouville equations 
ffTU|) for the marginal distribution functions F{t) = (^I , Fi(t) Fgit) we derive 

|^.(*) - ( - Efe. I") + ± - + 

1=1 ^ *^j=i 

^ f d d 

+ Y,J $(gi-g,+i),— )F,+i(t), (13) 

Fsit) \t=o= FsiO), s>l. (14) 



A solution of Cauchy problem f[T3l) - f[Ml) was constructed in [13|,[17|. The microscopic 
phase densities in this case are given by marginal microscopic phase densities ([1]) and their 
average values at time moment t G M 

oo 

{A){t) = {G{0),F{t)) =J2l dxi...dxsGs{0)Fs{t) = 
= [Git), F(0)) =Y,- dx,... dxsG,{t)Fs{0). (15) 
are determined by the expressions 

(G(^))(t) = F,(t,ei,...,es), s>i, 

where F^it, ^i, . . . , C,s) is the value of a solution of initial- value problem f[T5]) - ffT^ at a point 
^i, . . . ,^s- The evolution equations for the average values of microscopic phase densities will 
be considered in the next section. 



III. EVOLUTION OF AVERAGE VALUES OF MICROSCOPIC PHASE 
DENSITIES 

We derive the evolution equations of average values of microscopic phase densities and 
construct a solution of the initial- value problem of the obtained hierarchy of equations. 
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A. The BBGKY hierarchy for average values of the marginal microscopic 
phase densities 

We consider the general case, i.e. for the k-aij type microscopic phase density. According 
to (fT^ and from the dual BBGKY hierarchy we derive 

i=l i^j=l 

+ E / ^e.+i(^$(r.-r,+i),^)(G('=+^))(t), (16) 
1=1"' * * 

with the initial data 

(G('=))(t,ei,...,a)|t=o=(G(^))(0), k>l. (17) 
Due to functional ( fT5l) . initial data (ITTIl are given as the functions 

(G('^))(0,ei,...,efc)=i^fc(0,ei,---,efc), 

where -Fs(0, .^i, . . . , ^s) is the value of the initial marginal state at a point . . . , ^s- 

We note that, according to the definition of functional (fT2!) and (ITSl) . the equality 
{G^^^){t) = {N^''^){t) holds in the case of finitely many particles. In the thermodynamic, 
limit the value {N^''^)(t) tends to (G'(*''))(t), i.e. to the solution of Cauchy problem (fT6l)-(fT71). 



B. The BBGKY hierarchy for average values of the microscopic phase densities 



For comparison, we derive hierarchy f[T6|l from equations ffTOj) . In the general case, i.e. 
the /c-ary type microscopic phase density (jll), we have 

i=i i^j=i * 

+ E / ^e.+i(^$(r.-r,+0,^)(iV^'+'^)(^), (18) 
with the initial data {k > 1) 

oo „ 

(iVW)(t,ei,...,efc)k=0= (1,^(0))"'5^^ J dXk+i...dXk+nDk+n{0), (19) 



n=0 
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where Dk+n{0) = Dk+n{0, ^i, . . . , ^k', Xk+i, • • • , Xk+n) is the value of the distribution function 
Dk+niO) of the initial state at a point ^i, . . . , ^k] Xk+i, ■ ■ ■ , Xk+n- 

As a result of the formal transition to the thermodynamic limit, Cauchy problem (I18I) - (IT^ 
gets form ([TB])-([T71). 

Let us transform equations (1181) to the form which is usually used in the plasma theory. 
We find the covariation of the microscopic phase density {5N'^^\^i)5N^^\^2))^ where the 
observable ^iV^^) = A^^^) - (iV^^)) is a fluctuation of the microscopic phase density. 
Since (5Ar(i)(ei)5iV(i)(6)) = {N'^^\ii)N^^\i2)) - {N^^\ii)){N'^^\i2)) , we have 

(Ar«(ei)ArW(e2))(t) = 

oo ^ „ n n 

= (1, D{Q)y^ 5Z / • • • dXnDn{Q, X^,...,Xn)Y,Y. ^(^1 ' MtM^2 " X,{t)). 

n=0 ' 1=1 j=l 

In accordance with the equality 

n n n n 

'^(^1 - ^^)^(^2 - Xj) = ^i^l - -X,)+J2 '^(^1 - ^O'^fe - Xj), 

i=l j=l 1=1 iy^j=l 

we find 

Thus, the covariation of the microscopic phase density {5N^^\^i)5N^^\^2)) is defined as 
follows: 

(5Ar«(ei)5Ar«(6))(t) = 5(6 - 6)(iV^'))(t,6) + (iV^'^)(t,ei,6) - (Ar«)(t, 6)(Ar(i))(t, 6)- 

For the regularized interaction potential, i.e. $'(0) = 0, in terms of the covariation of 
the microscopic phase density {5N^-^\^i)5N^^\^2)) , the first equation from hierarchy f|T8|) 
reduces to the Vlasov-type equation 

|(iv(^^)(t,ei) + (^^i,^)(iv(^))(t,ei)- 

= / ^6(^$(ri -r2),^)(5Ar«(6)5Ar«(6))(t)- (20) 

This unclosed equation with respect to the average value of additive-type microscopic phase 
density {N^^^)(t) is the conventional equation for the description of a system of charged 
particles 
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C. The system of charged particles 



If we consider a plasma, i.e. a system of charged particles, then $(ri — = |ri — ^, 
and the macroscopic electric field {E){t) is defined from the equation 

div(E)(t,0 = e(iV«)(t,0 = 

oo 

where Nn\t,C,) is the additive- type microscopic phase density (jll). 
Thus, in this case, equation fl20l) takes the form 

Here 6E = E — (E) is a fluctuation of the electric field E = {Ei, . . . , En, . . .), where 

_ p., ^ d f Ni'\t,^';x,,...,Xn) 

En = E{t, r;xi,..., Xn) = -e tt / 



and {E) = {E){t,r) is the macroscopic electric field 



J I ■y^ I 

Equation (12T1) is an unclosed equation with respect to the average value of the microscopic 
phase density {N^^^)(t). 



D. The hierarchy for average values of the microscopic phase densities: 
the evolution of states 

According to f|T2|) and Cauchy problem f|TO|) - f|TT]) of the Liouville equations the evolution 
equations of the average values of the k-aij type microscopic phase densities get the form 
of hierarchy f|T8l) with initial data f|T9l) . 

Using the BBGKY hierarchy ( fT3l) . due to formula (fT5l) for marginal phase densities ([T]) 
we derive hierarchy ( |T6i) with initial data ( |T7I) . 

As we noted above, hierarchy (JT8l) in the thermodynamic limit transforms to the hierarchy 
of equations (fT6|) . 
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E. On solution of the initial-value problem 



To determine a solution of hierarchy f|T6l) . we introduce some prehminaries. 
On integrable functions, we define the following group of operators 

{Sk{-t)fk){Ci,---,Ck) :=/,.(Si(-t,ei,...,efc),---,Sfc(-t,ei,...,6)), (22) 

where the functions = (Vi(t), are the solution of the Cauchy problem of the 

Hamilton equations for "macroscopic variables" (dynamics of continuum) 

|fl.(<) = V.(t}, 

s^-.w = - t gi^f(m-m)). (23) 

with the initial data 

Ri{0) = n, 

Vi(0) =Vi, i = l,...,k. 

The generator of group fl22l) is defined by the Poisson bracket with respect to the variables 
^1, ■ ■ ■ on the continuously differentiable functions fk = fki^i, ■ ■ ■ , ^k) 

lim i (S,(-t) -/)/.= (- h-., l) + t i^^i-. --■.). |-» A- 

A solution of Cauchy problem (|T6l) -( fT7ll is defined by the expansion over the arity index 
of the microscopic phase density, whose evolution is governed by the corresponding-order 
cumulant (semiinvariant) of the evolution operators fl22l) . namely 

(G('=))(t,^i, . . . ,e.) = (f/(t)(G)(0))^(ei, . . . = 

oo „ 

= ^^. / ^efc+l---rfefc+n2tl+„(-t,ri,6+l,...,efc+n)(G'('=+"^)(0), (24) 
n=0 ' 

where (G')(0) = (0, {G^^^){0), . . . , (G'(''))(0), . . .) is the sequence of integrable functions. If 

n > 0, 

2ti+„(-t,yi,efc+i,...,efc+n)= 5^ {-lf^-\\P\-l)\l[ S\x,\{-t,X,) (25) 

P:{yi,x\y}=u,x, x,cP 
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is the (1 + n)t/i-order cumulant of the groups of operators (122|) . is the sum over all 
possible partitions P of the set {Yi,C^k+i, ■ ■ ■ ,^k+n} into |P| nonempty mutually disjoint 
subsets Xi C {Yi,X \ Y} = {Yi,C,k+i, ■ ■ ■ ,^k+n}- The set Yi consists of one element of 
Y = . . . , ^fc), i.e. the set (^i, . . . is a connected subset of the partition P (|P| = 1). 
If (G*-^'')(0) are integrable functions |17|], series fl24p converges for small densities (see also 

Q,Q)- 

We remark that a solution of Cauchy problem f|T6l) - f|T7j) can be represented as the iteration 
series of the BBGKY hierarchy ( [T6|) 

(G'(^))(t, ei, . . . , a) = E / • • • / / ■ ■ ■ ^^'^+" ^ 

n=o{ { 

k 

xSk{-t + tl) 5^ ( - Ant(^l, k + l))Sk+l{-h +t2)... 

k+n-1 

...Sk+n-l{-tn-l+tn) J] ( - Ant (^n, ^ + n)) (G^'-""^ ) (0) , (26) 



«n = l 



where 



-£>„.(i.j)-(|-4(r<-r,),|->. 

If we apply the Duhamel formula to cumulants of groups of operators fl2^ . solution expansion 
reduces to iteration series fl26p of the BBGKY hierarchy. 



IV. EVOLUTION EQUATIONS OF NON-MARKOVIAN TRANSPORT 

Using above-obtained results, we develop a new approach to the description of the evo- 
lution of the average value of microscopic phase density. Namely we formulate the non- 
Markovian evolution equation for the average value of microscopic phase density. 



A. The generalized equation for an average value of the microscopic phase density 

Consider Cauchy problem f|T6l) - f|T7j) with initial data which completely defined by the 
average value of additive- type microscopic phase density {G^^^){0), for example, 

k 

(G('=))(0) = n(G'(^^)(0,e.). (27) 

i=l 



13 



We note that {G^^^){0, ^i) = Fi{0, ^i). Therefore, initial data (1271) have the transparent sense: 
they satisfy the chaos condition. 

In that case, the initial- value problem of the BBGKY hierarchy (fT6|) - (fT7|) is not a 
completely well-defined Cauchy problem, because the generic initial data are not inde- 
pendent for every unknown functions {G^^'^){t), A; > 1, of the hierarchy of equations. 
Thus, it naturally arises the possibility of reformulating such initial-value problem as a 
new Cauchy problem for independent unknown function, i.e. the average value of the 
additive-type microscopic phase density {G^^^){t), together with explicitly defined function- 
als (G'('=))(t,^i,...,^fe I k>2, of the solution {G^^^){t) of this Cauchy problem 

instead other unknown average values of microscopic phase densities [l^, jl8 |. 

The functional (G'W)(t,G, •••,61 k > 2, are represented by the expansions 

over the products with respect to the average value of the additive-type microscopic phase 
density {G^^^){t) 

OD ^ „ k+n 

{G^''^){t,Ci,...,^k I (G(i^)W) ■■=J2^_ / d^k+i. . .dCk+n'Si+n{t)l[{G^'^){t,^i), (28) 

n=0 ' i=l 

where the evolution operators 5Ji+„(t) = 5Ji4.„(t, Yi, 6+i, . . . , 6+„) are defined from the 
condition that functional {G^^''){t)[t \ {G^^^){t)) are congruent with solutions (|24l) of initial- 
value problem (fT6|) - (fT7|) . 

We give examples of first terms of series (|28|) : 



QJi(t,yi) = 2ti(t,yi 



where Yi = (6, •••,6) and 2tfc(t) is the kth-oidei, k > 1, cumulant of the scattering 
operators 

k 

Skit) = Sk{-t,Ci, • • .,^k)l[Si{t,^,), (29) 

1=1 

the operator Sk{—t,^i, . . . ,$,k) is defined by formula ( l22ll . and Si(t) = I is the identity 
operator. 

On integrable functions, the action of scattering operators (1291) is defined by the formula 
{Skit)fk) (6, • • • , a) = fk (Si {t, Si(-t, ei, • . . , a)) , • • • , Sfc (t, Eki-t, 6, • • • , 6))) , (30) 
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where the functions Sj(t) = {Vi{t), i?j(t)), i = 1, . . . ,n are solutions of the Cauchy problem 
of Hamilton equations (123|) for "macroscopic variables" with corresponding initial data. The 
generator of group fl29|) of scattering operators is defined by the Poisson bracket with respect 
to the variables ,^1, . . . , with an interaction potential $ on the continuously differentiable 
functions fk = fk{^i,..., ^k) 

iinji(s.(t)-/)A=x:(|-*(''.-o). 

In terms of scattering operators (l29l) . the first terms of series (l28l) have the form 
^^{t,Y,)=Sk{t,Y), 

k 

QJ2(t,l^i,a+i) = - Sk{t,Y)Y.^2{t,^j,^k+i) + {k- l)Sk{t,Y). 

i=i 

If {G^^^){t) is an integrable function, then, in case of low densities, expansion (!28|) is the 
converging series [isl l. 

The average value of the additive- type microscopic phase density {G''^^){t) is governed by 
the following Cauchy problem {the generalized evolution equation for the average value of 
microscopic phase density): 

= I de2(^<f(ri -r2),^)(G(2))(t,ei,6 I (31) 
with the initial data 

(G«)(t,6)|t=o = (G«)(0,6)- (32) 

The functional (G*^^-') (t, ^1, ^2 | {G^^^){t)) in the collision integral of equation (l3Ti) is defined 
by expansion (1281) 

2 

(G(^^)(t)(t,ei,6 1 (G«)(t)) =2ti(t,ri)n(G'(^^)('^,e.) + 

i=l 

/s 3 
f/e3(2i2(t,Fi,a+i)-2ti(t,Fi)^2l2(t,0,a+i))n(^^'^)('^'^^) + ---' (33) 
j=i i=i 

where Fi = (^1,^2)- 
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We represent the first term of expansion (133|) in an explicit form. According to the 
definition of two-particle scattering operator (l30l) . we have 

2 

%^{t,Y^)l[{G^''>){t,^,)= (34) 

i=l 

= Si(t, Si(-t, ei,6)))(G'«)(^' 6,6))), 

where the functions (i = 1, 2) 

Ei{-t, 6,6)) = 6,6),^.(-i, ^1,^2) +tVi{-t, 6,6)) 

are defined as above by (123!) in formula (!30|) . 

The following statement is true. If the initial data are completely defined by (G''-^^)(0), 
then Cauchy problem (fT6ll - (fT7l) is equivalent to initial- value problem (j3Tl) -( !32l) of the gen- 
eralized evolution equation for the average value of the microscopic phase density and the 
functionals {G'^''^) (t, 6, • • • , 6 I {G''^^){t)), k >2, defined by expansions provided a low 
density of particles. 

Thus, at an arbitrary moment time the evolution of the average value of the microscopic 
phase density can be described by equation (!3T!) without any approximations. 
The solution of Cauchy problem ( l3Tl) -( !32l) is defined by the expansion 

00 „ 71+1 

(G('^)(i,6) = E;^ / rf6...c?6+n2li+„(-t,6,...,6+i)n<^^'^)(0'^^)' (35) 

n=0 ■ i=l 

where 2ti+„(— t, 6, • • • , 6+1) is the (1 + n)th-oTdeT cumulant of the group of operators (!22ll . 

For the low densities, i.e. J \{G^^^){0,$,)\d^ < e~^, series ( l35l) converges. If {G^^^){0) is a 
continuously differentiable function with compact support, then, for t G M, it is a classical 
solution of generalized equation (!3T|) for the average value of the microscopic phase density. 

We observe the links of the introduced generalized equation for the average value of the 
microscopic phase density with Bogolyubov's method of the derivation of kinetic equations 
[igI i . Functionals (l28l) are formally similar to the corresponding functionals of Bogolyubov's 



method 191] if they satisfy the principle of weakening of correlations [13]. The proof of this 
statement is completely similar to the proof of the equivalence of both representations of 
the BBGKY hierarchy solutions as iteration series (|26|) and functional series (l24l) . 
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B. An example: the regularized interaction potential 



We separate the Vlasov term in equation (!3T]) to represent it in the conventional form for 
the description of a system of charged particles. 

We introduce the macroscopic characteristic of correlations in continuum 

{G){t,^i,^2 I = I (G'«)W) - (G«)(t,6)(G'«)(t,e2), (36) 

where the functional (G'(2))(t) (t, ^i, 6 I is defined by expansion fl33p . Then the 

generalized evolution equation (^11) for the average value of microscopic phase density gets 
the form 



+ / c^6(^$(n -r2),^)(^)(t,ei,6 I {G^'^m) (37) 
with initial data ( l32i) . 

In the case of regularized interaction potential, i.e. $'(0) = 0, the functional (Q) (t, ^i, ,^2 I 
{G^^^){t)) is the formal thermodynamic limit of the covariation of the microscopic phase 
density {5N^^\^i)5N^^\^2)) from the right hand side of equation fl20l) . 

If we consider a plasma, i.e. a system of charged particles, the macroscopic electric field 
(E) (t) is determined from the equation 

div(E)(t,0 = e(G«)(t,0. 
In this case, equation (1371) gets the form 

(| + -"^ + »(£^")(*.'-)^^)(G'")(«) = / J^jT^Ti a|;(e>(«.«'l <G'")W). 

where the functional {G)[t \ {G^^^){t)) in the collision integral is defined by expansions 



C. On the Fokker-Planck representation of a generalized collision integral 



If we consider a plasma IGj, then the small parameter is the plasma parameter associated 
with the weak interaction 
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where r^) is the Debye screening parameter, e is the electric charge and f3~^ is the value of 
the order of the average kinetic energy of electrons. 

Let us consider a first term of the collision integral expansion of (^)(t,^i,^2 I {G^^^){t)) 
in equation (^7}i . i.e. the first order term with respect to the plasma parameter associated 
with the weak interaction, 



(38) 



We will construct a suitable approximation of the non-Markovian collision integral fl38|) 
for a plasma. Due to the fact that, for the scattering operator S'2(t, 1,2), the Duhamel 
formula is formally true (see also flM|) ) 

t 

drS2H + r, 6, 6) ((^"^(^1 - ^2), ^) + 



+ (— <l>(ri - r,), j-))S,{-r + t, 6)5i(-r + t, 6), 
or2 0V2 

where 5*2 (—t + r, ^1,^2) is evolution operator (12^ . collision integral (1551) can be given as 
follows: 

t 



+ (^$(n - r2), + t,^^)S^i-T + t, 6)(G'«)(t, ei)(G'«)(t, 6). 

Then, by expanding the obtained collision integral X in the Dyson— Phillips series in 
the plasma parameter and by restricting ourselves to the first term of the expansion (weak 
coupling limit), we find 

t 



-r2 (-t + r)) , ( A _ ^) ) (t, ei) (t, 6) , 

where rj(— t + r) = 'rj(— t + t, ^1, ^2), = 1, 2, are the solutions of initial-value problem (123|) 
for k = 2. We remark that this expression has similar structure to the generalized Landau 
collision integral in the kinetic theory [I61] . 
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Thus, we can finally rewrite equation (137|) in the Fokker-Planck representation [20| as 



where the transport coefficients are defined by the expressions 

t 



p-/^(t) = j dr j d^2^$(ri-r2)^<l>(ri(-t + r)-r2(-t + r))(GW)(t,6), (40) 
^ ^ 

S"(t) = I rfr I rf6^$(ri - r^) (|^$(ri(-t + r) - r,[-t + r))^ - 



- A A<^'(n(-t + r) - r2(-t + r))) (G«)(t, 6)- (41) 
af or I ^ 

Evolution equation (!39|l is a non-Markovian evolution equation for the average value 
of the microscopic phase density with the Fokker-Planck collision integral determining by 
nonlinear transport coefficients (l40l) . (14T]) . It can be reduced to the canonical Fokker-Planck 
equation as a result of further approximations [3] ■ 

We observe that the higher-order terms in ( |3711 with respect to the plasma parameter from 
the functional {G){t, ^1,^2 \ {G^^^){t)^ are equal to zero in the weak coupling approximation. 

V. CONCLUSION 

The BBGKY hierarchy of equations f|T6l) for the generalized microscopic phase densities is 
formulated, and the microscopic derivation of a non-Markovian collision term in the evolution 
equation for the average value of the additive-type microscopic phase density ( !3T|) is done. 

The time-nonlocal generalization of the Fokker-Planck equation (l39l) -( HTl) on the basis 
of such collision term ( l38l) is proposed. The evolution equation generated by such collision 
term can be applied for further analysis of turbulent transport |7|. The memory effects can 



be important for the description of transport under saturated turbulence jgl. [lol. Ill|. 
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